Velocity of sound in relativistic heavy-ion collisions 
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We have studied the rapidity distribution of secondary hadrons produced in nucleus-nucleus col- 
lisions at ultra-relativistic energies within the ambit of the Landau's hydrodynamical model. A 
reasonable description of the data can also be obtained by using the Bjorken's hydrodynamical 
model if the boost invariance is restricted to a finite rapidity range. The sensitivity of the hadronic 
spectra on the equation of state vis- a -vis the velocity of sound has been discussed. The correlation 
between the velocity of sound and the freeze-out temperature has been indicated. The effects of 
the non-zero widths of various mesonic and baryonic degrees of freedom up to the mass value ~ 2.5 
GeV is seen to be small. 

PACS numbers: 25.75.-q,25.75.Dw,12.38.Mh 



I. INTRODUCTION 



One of the important problem in the field of rel- 
ativistic heavy-ion collisions is to find out the equa- 
tion of state (EOS) of the matter formed after nu- 
clear collisions at ultra-relativistic energies 0. Un- 
der the assumption of local thermal equilibrium, the 
EOS is the functional relation between pressure (P) 
and the energy density (e), where P and e are re- 
lated through the velocity of sound, Cs which is de- 
fined as = {dP/de)isentropic Q - For a massless, 
non-interacting gas, — 1/3 (ideal gas limit). The 
velocity of sound plays a crucial role in the hydro- 
dynamical evolution of the matter created in heavy- 
ion collision and affects, among others, the momen- 
tum distribution of the particles originating from the 
fiuid elements at the freeze-out stage. In the present 
work we will assume the order of the phase transition 
from QGP to hadrons to be first order. However, 
it should be mentioned here that this issue is not 
fully settled yet. The phase transition from QGP to 
hadrons could be weak first order, second order or 
it may be just a cross over depending on the mass 
of the dynamical quarks Consider a situation 
where quark-gluon plasma (QGP) is formed at the 
initial state. In such a scenario, the matter evolves 
from an initial QGP state to the hadronic phase via 
an intermediate mixed phase of QGP and hadrons 
due to the expansion of the system (hence cooling) 
in a first order phase transition scenario. Finally 
the system disassembles to hadrons (mainly pions) 
at the freeze-out where the interaction among the 
particles become too weak to maintain the equilib- 
rium. The velocity of sound is very different in the 
three stages of expansion mentioned above reflecting 
the interaction among the constituents of the mat- 
ter in the three stages. While in the QGP phase, it 
should in principle approach the ideal gas limit, in 
the mixed phase it should reduce to zero due to van- 
ishing pressure gradient, indicating "softness" of the 
EOS. Then below the critical temperature, it should 



have a value that reflects the presence of interacting 
hadrons in the system. 



Relativistic hydro-dynamical models have been 
routinely used to describe the multiplicity distribu- 
tion in the rapidity space, transverse mass spectra 
of hadrons etc produced in nuclear collisions. One 
of the aim of the present work is to calculate the ve- 
locity of sound for hadronic system and compare the 
results with lattice QCD calculations. In the present 
calculation we have taken into account the non-zero 
width of the unstable hadrons. We study the rapid- 
ity distribution of the particles within the framework 
of relativistic hydro-dynamical models proposed hy 
Landau and collaborators and by Bjorken |^. 
The results of the analysis performed within these 
models will be compared with experiments to extract 
the velocity of sound. It is expected that at AGS and 
also at SPS energies, the Landau hydro-dynamical 
model can be applied although the Bjorken hydro- 
dynamical model has been used at SPS energies. 
The main criticisms of Landau model are: (i) ne- 
glect of leading particle effects and (ii) removal of 
radiation energy due to the deceleration required in 
this model for full stopping. These difficulties, how- 
ever, can be removed if one assume that during the 
collisions the valence quarks move without much in- 
teraction and the energy carried by the gluon fields 
is stopped in the collision volume 6]. This assump- 
tion is justified because the gluon-gluon interaction 
cross-section is larger than quark-quark interaction 
due to larger color degeneracy of the gluons. The 
gluon field thermalizes after a time Ti, providing the 
initial condition of the Landau model. In this picture 
the removal of energy of the decelerated gluons fields 
due to the bremsstrahlung is prohibited by the color 
confinement mechanism. Under these conditions, it 
is obvious that only a fraction of the beam energy 
is stopped in the collisions, which can be taken into 
account by introducing an in-elasticity factor in the 
model iQj. 
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TABLE I: Particles taken for calculation. 



Baryons 




Mesons 








n 




r,(547-1440), rj' 


N(1440 - 


2600) /o(800-1710),/i(1285, 1420),/2 (1270-2340) 


A(1232 - 


2420) 


p(770-1700),p3 


A(1115 - 


2350) 


w(782-1650),L^3 






ao (980- 1450), ai 


E(1382 - 


1820) 


(^,<?!>3,7r(1300-1800),7r2 






/4,/ii ,7^1 (1270-1400) 


S(1530 - 


1820) 


a2,a4,bi,K2 









Bjorken hydrodynamical model predicts a plateau 
structure for the rapidity distribution of the fluid el- 
ements which is not observed experimentally for the 
entire rapidity range. It has been indicated in 
from perturbative QCD that the initial energy dis- 
tribution is a broad Gaussian in rapidity even at 
LHC energies. However, we will show that a good 
description of the data is possible at AGS and SPS 
energies in the framework of the Bjorken's model if 
we treat the upper limit of the fluid rapidity as a 
parameter. In other words the boost invariance has 
to be limited in a finite rapidity range for the de- 
scription of the data. 

The paper is organized as follows. In the next 
section we obtain the velocity of sound for a hadronic 
model and compare the result with the lattice QCD 
calculations |3] and those obtained from a simple 
confinement model as discussed in Ref . In section 
III we extract the velocity of sound from the particle 
number density in rapidity space at AGS and SPS 
energies. Finally in section IV we present summary 
and discussions. 



II. EQUATION OF STATE FOR HADRON 
GAS 

For the study of the EOS for the hadronic gas we 
take all the hadrons as listed in the particle data 
book up to the strange sector. The complete list 
is given in Table I. The thermodynamical quanti- 
ties like energy density (e), and pressure (P) can be 
calculated using the standard relations, 



g 



hadrons 



E{p)f{E)p{M) dlVr (Pp 



hadrons 



-9 

9 I P 



(27r)3 J ZE 



f{p)p{M)dM^d^p (1) 



where E — \/p^ + is the energy of the particle 
of three momentum p and invariant mass M, g is 
the internal degrees of freedom and f{E) is the well- 
known thermal distribution for bosons and fermions 



which is given by f{p) = exp{{E — p)/T) ± 1 

where /i is the chemical potential, T is the tempera- 
ture. The summation in Eq. ^ is carried out for all 
the hadrons up to strange sector (9j. The sensitiv- 
ity of the EOS on the hadronic and electromagnetic 
spectra has been studied in p^ . with less number 
of hadrons and the width of the hadrons are ignored. 
p{M) denotes the spectral function of the hadrons 
of pole mass m and width F given by. 



p{M) - - 



AfF(Af) 



TT (M2 - m?Y + M2F2 



(2) 



Note that for stable particle p{M) 6{M'^ - vn?) 
as F 0. To evaluate the spectral function of a 
hadron in a thermal bath one should consider all 
the elastic and in-elastic processes through which 
the hadrons under consideration interact with all 
the constituents of the thermal bath These 
interactions give rise to the momentum dependent 
effective mass and widths. However, in the present 
work we restrict to the vacuum values of these quan- 
tities which are taken from particle data book j^. 
The spectral function can also be related to the 
two body phase shift {5) of the various processes as 
p{M) cx d5{M)/dM, leading to the results obtained 
by S-matrix formalism of the statistical mechanics 
proposed by Dashen et al. (12| . 

From Eq. one can calculate the entropy density 
(s) by using the relation, s = (e + P — pn)/T where 
n is the net baryon density given by. 



E 

baryons 
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(27r)3 



f{E)p{M) dM^ d^p 



gef f is the effective degeneracy which is parameter- 
ized as geff{p,T) = 90s/(47r2T3). The velocity of 
sound for = can be calculated from the geff as. 



T ds 

7dr 



T dgeff 



9eff 



dT 



(3) 



We will assume the chemical potential to be zero 
for mesons and consider two cases for baryonic chem- 
ical potential, p, = and 200 MeV. 

In Fig^we plot the variation of various thermody- 
namic quantities: e , P, geff and as a function of 
temperature. We make two observations: (a) the ef- 
fects of the finite width in the masses of the hadrons 
on the thermodynamical quantities turn out to be 
small and (b) the effect of the finite baryonic chemi- 
cal potential leads to a slightly larger values of e , P 
and gef / but lower values of at (typical) freeze-out 
temperature of 120 MeV. 

For comparison of our results with those from the 
lattice QCD calculations Q, we parameterize the 
variation of energy density (e) of lattice QCD results 
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FIG. 1: Variation of energy density (e), pressure (P), ef- 
fective statistical degeneracy (geff) and velocity of sound 
(cs) with temperature (T). 



with temperature as follows Ts']: 



FIG. 2: Variation of energy density, g^ff and velocity of 
sound using the combined results of hadron gas, confine- 
ment model and lattice calculations with temperature. 
The lattice results for energy density has been read from 
01 . geff and Cs are derived by using the equations men- 
tioned in the text. 



e = T A tanh 



B{ — f 



(4) 



where A, B and C are parameters whose values are 
12.44, 0.517 and 10.04 respectively. Note that the 
effect of baryons in the EOS is neglected here. 

The evolution of the system under the assumption 
of boost-invariance along the longitudinal direction 
is governed by the equation 5], 



P 



0:P 



(5) 



where Cs is the velocity of the sound in the medium. 
We would like to mention here that the velocity 
of sound sets the expansion time scale {{Tf,xp)~^ ~ 
{l/e)de/dT = (1 + c^)/t) for the system. This time 
scale should be larger than the collision time scale 
i(Tcoii)~^ ~ nav, where a is the cross section, v is 
velocity and n is density) for thermal equilibrium to 
be maintained in the system. Therefore, determina- 
tion of velocity of sound becomes very important for 



the study of the space time evolution of the system 
in general and hadronic spectra in particular. 

We obtain the expression for the velocity of sound 
as, (see also il^), 



T 



1 



T/ cosh{B{^)C) sinh{B{^)C) 



(6) 

The results for the confinement model were ob- 
tained from by extracting QeffiT) by using the 
relation e/T^ = (7rV30)5e// and Eq.|Sl 

In Fig[21the energy density, effective statistical de- 
generacy and the velocity of sound evaluated in the 
hadronic model at zero baryonic chemical potential 
are compared with the lattice data (/i — 0) and those 
obtained form the confinement model. We observe 
that, (a) our results match with those from the lat- 
tice and the confinement model near Tc, (b) they are 
higher than those obtained from lattice below Tc (In 
this context we would like to point out two things. 
Firstly, lattice results usually employ quark masses 
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which are too large, leading to larger pion masses, 
^lat > ^^physicai ^ r^^iQ resulting thermal suppres- 
sion of these degree of freedom causes a considerable 
discrepancy for lattice EOS with respect to EOS for 
hadron gas. Secondly, the lattice calculations below 
Tc has large errors), (c) at temperatures below Tc 
the velocity of sound shows an interesting trend, it 
increases with decrease in temperature then falls to 
zero as temperature of the system approaches zero 
and (d) the complete EOS for the system with ini- 
tial QGP state, may be those obtained from lattice 
for temperatures above Tc and those given by the 
present calculations for temperatures below Tc- 

The velocity of sound, which is an input to the 
hydro-dynamical evolution of the system via EOS 
influence the rapidity distribution of the hadrons. 
In the next section we try to extract the velocity 
of sound from the rapidity distributions of various 
hadrons produced from heavy-ion collisions for vari- 
ous collision energies using hydro-dynamical model. 
We then compare the values obtained with the re- 
sults of the model under consideration. 



III. RAPIDITY DISTRIBUTION OF 
SECONDARIES AT AGS AND SPS 
ENERGIES 



TABLE 11: AGS results 



Type 



X 



Probability 



proton 4 AGeV 1/5 (1/3) 0.65 (34.5) 1.0 (0.023) 

proton 6 AGeV 1/5 (1/3) 2.15 (40.2) 1.0 (0.010) 

proton 8 AGeV 1/5 (1/3) 7.1 (74.8) 0.999 (0.00001) 

proton 11.6 AGeV 1/5 (1/3) 1.49 (31.1) 1.0 (0.0388) 



model is given by 0, Q , 

^ = -nRHsof3csexp[f3u;f] \lo{q) - ^h{q) (7) 
ay L q 



where, q — ^^uP'^ — c^y"^, ujf = ln{Tf /Tq), Tf is the 

freeze-out temperature, y is the rapidity, R is the 
radius of the nuclei, 21 is the initial length, sq is the 
initial entropy density, 2(3 — {1 — cj.)/c^ and Iq, Ji 
are the Bessel's function. The quantity ttR^Isq is 
fixed to normalize the experimental data at the mid 
rapidity. For ujf » Cgy the quantity dS/dy can be 
approximated by a Gaussian distribution. 



dS exp(- 

— ^ Const. -jzz 

dy V27r(T- 



> 



(8) 



where a = 2cj//(l — c^). 

The entropy density in the Brojken hydrodynam- 
ics is given by [31 , 



There are two well known hydro-dynamical mod- 
els as mentioned before 0, Q for the description of 
the space time evolution of the system formed af- 
ter the collision of heavy-ions. It may be mentioned 
that in the coherent interactions the collective effect 
is an important feature, unlike incoherent collisions 
where collision is considered to be a succession of 
independent nucleon-nucleon interactions. The ba- 
sic difference is, in Landau's model, in the center of 
mass of the collision two nuclei would be stopped 
and all the kinetic energy would be used up for the 
particle production, while the baryon transparency 
in the mid-rapidity region is the basic feature of the 
Bjorken's model. However, the Landau's model can 
be applied to describe the system formed after nu- 
clear collision with appropriate initial conditions as 
mentioned in the introduction. It is expected that 
the Landau hydrodynamics will work well at AGS 
and may be at SPS energies, while the Brojken hy- 
drodynamics should in principle work well for ener- 
gies at RHIC and LHC. 

In the following we briefly discuss the rapidity 
distribution of the particles produced in relativis- 
tic heavy-ion collisions in Landau's and Brojken's 
model. 

The amount of entropy {dS) contained within a 
(fluid) rapidity dy in the Landau hydro-dynamical 



dS_ 
dy 



Tf 



Sfexp[—2l3ijjf\ 



(9) 



where, A is a constant and s/ is the entropy density 
at the freeze-out. It is to be noted that the quantity 
dS/dy is independent of the rapidity in accordance 
with the assumption of boost invariance along the 
longitudinal direction 'H'l . 

The multiplicity distribution of the secondaries is 
obtained by folding the multiplicity density in rapid- 
ity space mentioned above by the thermal distribu- 
tion of the fluid elements, 



dN 
dY 



(27r)3 J dy 



dN 



/{Mt, y, Y)MTCOsh{Y - y)dyd^pT 

(10) 

where, dN is the number of particles within the 
rapidity interval dY , Mt is the transverse mass 
(= ^Jp\ + A'P) of the particle, E = MTCosh{Y-y), 
is the energy and f{E) = l/{exp[E/T] ± 1). dS/dy 
is related to dN/dy by a constant factor (see for 
details). Performing the px integration in the above 
equation we obtain, 

^ « y -^Ky^ Y'l ™7 T) exp{-mcosh{Y-y)/T) dy 
^ (11) 
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FIG. 3: Rapidity spectra for protons at 4,6,8, and 11.6 
AGeV Au+Au collisions compared to rapidity spectra 
obtained from Landau hydrodynamics with velocity of 
sound 0.2 (solid line) and 0.333 (dashed line). 
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FIG. 4: Rapidity spectra for pions at 40,80,158 A GeV 
Pb-|-Pb collisions compared to rapidity spectra obtained 
from Landau hydrodynamics with velocity of sound 0.2 
(solid line) and 0.333 (dashed line). 



h{y,Y;m,T) — — ^+2 — -H —r- - 

T coshiY — y) cosh^yY^y) 

(12) 

It may noted from Eqs^] and E| that the parti- 
cle mass tend to make the rapidity distribution nar- 
rower. 

It may be mentioned that, while in Landau's 
model the width of the rapidity distribution is sensi- 
tive to the velocity of sound and the freeze-out tem- 
perature, in the case of Brojken's model it is not (be- 
cause of its independence on rapidity, the quantity, 
dS/dy can be fixed here by the normalization at the 
mid-rapidity). While integrating over y in Eq. IIII 
we treat the range of ?/ as a parameter in case of 
Bjorken's hydrodynamics In case of Landau's 
hydrodynamics the integration limit for rapidity is 
infinite. We apply the Landau's model to the rapid- 
ity distributions of the produced hadrons at AGS 
and SPS energies. It will be seen later that the Lan- 
dau model fits the pion rapidity spectra well at lower 
SPS energies, however at highest SPS energy the de- 
scription is not very satisfactory, therefore we do not 
attempt to apply the model to RHIC energies. We 
fix the freeze-out temperature Tf to be 120 MeV a 
value obtained by studying the transverse momen- 
tum distributions of the hadrons 0, 0| . 

The rapidity distribution of the protons 0, 
at AGS energies are well described by the Lan- 
dau hydro-dynamical model with velocity of sound, 

= 1/5, a value different from that corresponding 
to an ideal gas (FiglS)). For — 1/5 the values of the 



TABLE III: SPS results for pions 
Type Probability 
TT" 40 AGeV 1/5 (1/3) 1.03 (51.7) 1.0 (0.0013) 
TT- 80 AGeV 1/5 (1/3) 2.6 (66.3) 0.999 (7 x 10"^) 
■K- 158 AGeV 1/5 (1/3) 19.0 (146.0) 0.80(4 x 10'^'^) 

chi-square is considerably smaller than for = 1/3 
(Table Inj) . The peak of the distribution at the mid- 
rapidity in the experimental data seems to indicate a 
large deposition of collision energy in these interac- 
tions, where the Landau's hydro-dynamical picture 
is applicable. 

For SPS energies the density of pions in rapidity 
space I23 is well reproduced (Fig^ within the am- 
bit of Landau's model with = 1/5 and Tf = 120 
MeV for 40 and 80 GeV/ A beam energies. For 
Ebea-ni = 158 GcV/A the worsens (Table HID- 
The Bjorken's model also does not give a satisfac- 



TABLE IV: SPS results for kaons 

Type (2 Probability 

K+ 40 AGeV 1/5 (1/3) 0.463 (18.0) 1.0 (0.38) 
K- 40 AGeV 1/5 (1/3.3) 0.56 (0.16) 1.0 (1.0) 
K+ 80 AGeV 1/4.5 (1/3) 0.4 (12.7) 1.0 (0.75) 
7^-80 AGeV 1/5 (1/3.3) 3.0 (0.67) 0.99 (1.0) 
K+ 158 AGeV 1/4.5 (1/3) 0.39 (15.6) 1.0 (0.55) 
K- 158 AGeV 1/5 (1/3.5) 3.8 (0.8) 0.99 (1.0) 
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FIG. 5: Rapidity spectra for pions at 40,80,158 A GeV 
Pb+Pb collisions compared to rapidity spectra obtained 
from Bjorken's hydrodynamics. The value of ymax{= 
—ymin) used iu the integration of Eg. 1111 are 1.44, 1.65 
and 1.6 for beam energies 40 AGeV, 80 AGeV and 158 
AGeV respectively. 



FIG. 6: Rapidity spectra for K+ at 40,80,158 A GeV 
Pb+Pb collisions compared to rapidity spectra obtained 
from Landau hydrodynamics with velocity of sound 0.2 
(solid line) and 0.333 (dashed line) and Tf = 120 MeV. 



tory description of the data at higher SPS energies 
(158 AGeV). However, as shown in Fig. [Sjat lower 
energies the Bjorken's model give a reasonable de- 
scription of the data. 

Interestingly, the spectra is well reproduced 
(Figs El and [7J) with the same freeze-out tempera- 
ture and velocity of sound as mentioned above for 
all the beam energies for both the models. How- 
ever, the spectra (Fig|SJ) is reproduced with a 
slightly higher freeze-out temperature, 132 MeV. A 
good description of the data for 40 and 80 GeV/ A 
beam energies is also obtained for T/=120 MeV and 
Cg = 1/3. 5(x^ = 0.1 and 0.3 respectively). For 
Ebeam = 158 GcV/A wc obtain cl = 1/4 for same 
freeze-out temperature (x^ = 0.3). The values of 
are shown in Table IIVI It is well-known that 
the system formed in nuclear collisions at AGS and 
SPS energies has non-zero baryonic density and in 
such a situation the differences in the values of the 
freeze-out parameters for and K~ are expected 
(see for details). 

The pion data from SPS at 40, 80 and 158 AGeV 
energies can also be reproduced if we use the Gaus- 
sian approximation of Eq. |S1 in Eq. 1111 and treat 
the width of the Gaussian (cr) as a parameter. Re- 
sults are shown in Fig|51 The value of cr = 0.9, 1.1 
and 1.25 for pion at beam energies 40, 80 and 158 
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FIG. 7: Rapidity spectra for kaons at 40,80,158 A GeV 
Pb+Pb collisions compared to rapidity spectra obtained 
from Bjorken's hydrodynamics. The value of j/maa;(= 
—y-min) used iu the integration of Ea. lTTl are 1.4, 1.5 and 
1.6 for beam energies 40 AGeV, 80 AGeV and 158 AGeV 
respectively. 
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FIG. 8: Rapidity spectra for K~ at 40,80,158 A GeV 
Pb-|-Pb collisions compared to rapidity spectra obtained 
from Landau hydrodynamics with velocity of sound 0.2 
for T/ = 120 MeV (dashed line) and 132 MeV (solid 
line). 



FIG. 9: Rapidity spectra for pions at 40,80,158 A GeV 
Pb-fPb collisions compared to rapidity spectra obtained 
from EgslSlandllll The width of the Gaussian in Eq. |S1 
is a parameter here. 



AGeV respectively. Results obtained for kaons in 
this procedure is shown in Fig. 1101 It is important 
to note that the values of the widths of the Gaus- 
sian which represent the rapidity distribution of the 
fluid elements are 1.25, 1.45 and 1.95 for 40, 80 and 
158 AGeV energies respectively and these values of 
a are consistently higher than the values obtained 
for pions. The difference can be attributed to the 
mass difference between pions and kaons because pi- 
ons and kaons are subjected to the same longitudinal 
flow velocity. 

In Fig. we show the constant contour in the 
— Tf plane. The results indicate that we need a 
value of less than 1/3 for a reasonable descrip- 
tion of the data. We observe that ~ 1/5 near 
the freeze-out for different collision energies and for 
various hadronic species. This indicates a certain 
kind universality for the hadronic matter produced 
in heavy-ion collisions at the late stage of the evolu- 
tion (freeze-out). 



IV. SUMMARY 

We have evaluated the velocity of sound in a 
hadronic model with hadronic degrees of freedom up 
to strange sector. This calculated velocity of sound 
is compared with the results obtained from lattice 



QCD calculations j^] and those obtained for the con- 
finement model m. The values compare well at the 
critical temperature. The effects of non-zero width 
of hadrons on the EOS is found to be small. Tak- 
ing a non-zero baryonic chemical potential slightly 
decreases the velocity of sound. Fixing the value of 
the freeze-out temperature (^ 120 MeV) from the 
•pT spectra of hadrons [TgL [l7j | we find that a value 
of = 1/5 give a good description of the data for 
Landau's hydrodynamical model. The data is also 
well reproduced within the ambit of the hydrody- 
namical model proposed by Bjorken if the boost in- 
variance is restricted to finite rapidity range. The 
value of = 1/5 indicates that the expansion of 
the system is slower in comparison to the ideal gas 
scenario (c^ = 1/3) and hence for such a system the 
maintenance of thermal equilibrium becomes easier. 
We observe that different hadron species produced 
in the nuclear collisions at different energies are well 
reproduced by a value of ~ 1/5, indicating some 
kind of universality of the matter at the freeze-out 
stage. 

Several other effects which are ignored in the 
present work need to be mentioned at this junc- 
ture. For a complete description of the space time 
evolution of the system formed after heavy ion col- 
lisions both the transverse and longitudinal expan- 
sion should be considered. The transverse momen- 
tum spectra of the hadrons is strongly affected by 
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the transverse flow. The rapidity distribution of 
the hadrons which is obtained after integrating out 
the transverse momentum may not be strongly af- 
fected by the transverse flow. However, the normal- 
ization of the rapidity spectra may be substantially 
affected by the transverse flow. The other quan- 
tities e.g. which may affect the normalization are 
non-zero chemical potential of the mesons, arising 
from the lack of chemical equilibrium in the system. 
The hadrons (pions, kaons and protons) originating 
for the decays of various mesonic {p, uj, (j) etc) and 
baryonic (A etc) resonances are ignored here. If the 
distribution of the resonances is homogeneous in the 
rapidity then the effects of the hadrons originating 
from their decays on the rapidity distributions can 
be ignored, otherwise these effects may change the 
distribution depending on the degree inhomogeneity 
and the abundances of the resonances at the freeze- 
out point. With the centrality of the collisions the 
normalization and the width of the distribution of 
hadrons in rapidity space changes. The normaliza- 
tion in the present work is treated as a parameter 
and the change in the width with centrality is rather 
smaU. At 158 AGeV SPS energy the width of the ra- 
pidity distribution of charged particle increases from 
1.5 to 1.62 when the centrality changes from (0—5)% 
to (25- 35)% H. 
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